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ABSTRACT 



Context. Infrared dark clouds (IRDCs) are promising sites to study the earliest formation stages of stellar clusters and high-mass stars, 
and the physics of molecular-cloud formation and fragmentation. 

Aims. The aim of this study is to better understand the physical and chemical properties of the filamentary IRDC G304.74+01.32 
(hereafter, G304.74). In particular, we aim to investigate the kinematics and dynamical state of the cloud and clumps within it, and 
the amount of CO depletion. 

Methods. All the submillimetre peak positions in the cloud identified from our previous LABOCA 870-/jm map were observed in 
C'''0(2 - 1) with APEX. These are the first line observations along the whole filament made so far. Selected positions were also 
observed in the ''C0(2 - 1), SiO(5 - 4), and CHsOHCSt - Ak) transitions at ~ 1 mm wavelength. 

Results. The C"0 lines were detected towards all target positions at similar radial velocities. CO does not appear to be significantly 
depleted in the clumps; the largest depletion factors are only ~ 2. Two- to three methanol 5i,-4k lines near ~ 241.8 GHz were detected 
towards all selected target positions, whereas SiO(5 - 4) was seen in only one of these positions, namely SMM 3. In the band covering 
SiO(5 - 4), we also detected DCN(3 - 2) line towards SMM 3. The '^^C0(2 - 1) lines show blue asymmetric profiles, indicating large- 
scale infall motions. The clumps show trans- to supersonic non-thermal motions, and virial-parameter analysis suggests that most of 
them are gravitationally bound. The external pressure may also play a non-negligible role in the dynamics. The analysis suggests that 
the fragmentation of the filament into clumps is caused by "sausage"-type instability, in agreement with results from other IRDCs. 
Conclusions. The uniform C"0 radial velocities along the G304.74 cloud shows that it is a coherent filamentary structure. Although 
the clumps appear to be gravitationally bound, the ambient turbulent ram pressure may be an important factor in the cloud dynamics. 
This is qualitatively consistent with our earlier suggestion that the filament was formed by converging supersonic turbulent flows. The 
poloidal magnetic field could resist the radial cloud collapse, which conforms to the low infall velocites derived. The cloud may not 
be able to form high-mass stars on the basis of mass-size threshold, but "only" stellar clusters and/or intermediate-mass stars. The 
star-formation activity in the cloud, such as outflows, is likely responsible in releasing CO from the icy grain mantles back into the 
gas phase. Shocks related to outflows may have also injected CH3OH, SiO, and DCN into the gas-phase in SMM 3. 
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1. Introduction 
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Infrared dark clouds (IRDCs) are seen as dark absorption fea- 
tures against the Galactic background radiation, particularly at 
mid-infrared wavelengths (IPerault et al. 19961 |Eganetal. 1998 



> ■ ISimon et al. 20061 IPeretto & Fuller 20091) . 



X 



Since their dis 

covery in the mid to late nineties, and after the first detailed 
studies ( [Carey et al. 1998] 2000), IRD Cs have been the 
subject of great research interest (e.g., ITeyssier et al. 2002 
■Johnstone et al. 2003 ; Rathborne et al. 2005 , 2006; 



iPillai et al. 2006. .Beuther & Sridhara n 2007: Sakai et al. 20081 

[Jimenez-Serra et al . 20101 
IKainulainen et al. 201 lb. and many 



Ragan et al. 2009 



Leurini et al. 201 



more). One of the reasons for such an active study of IRDCs 
is that they are found to be promising targets to study the 
initial conditions and early stages of high-mass star (> 8 



M0) and stellar cluster formation, both of which are quite 
poorly known compared to the formation of solar-type stars. 
Moreover, IRDCs are interesting objects for the studies of 
molecular-cloud formation and their fragmentation into clumps 
and cores [e.g., [Miettinen & Harju 2010 (hereafter. Paper I); 
IJackson et al. 2010r Wang et al. 201 1| . Infrared dark clouds are 
often found to be filamentary in shape, resembling the mor- 
phologies seen in numerical studies of cloud formation (e.g., 
Nakajima & Hanawa 1996[ IBallesteros-Paredes et al. 19991 
Klessen & Burkert 20001 IKlessen 20011 IHeitsch et al. 20081 
2011). On the other hand, theories and models of insta- 
bility and fragmentation of filamentary, or cylindrical, 
gas structures have been extensively studied in the past 
(e.g., IChandrasekhar & Fermi 1953] IStodoUciewic z 19631 



lOstriker 19641 IBastien 19831 [Tnutsuka & Miyama 1992 
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Curry 2000). The filamentary IRDCs are well-suited for testing 
these theories/models, such as their predictions about the 
preferred fragmentation length-scales. 

The target cloud of the present study is the IRDC 
G304.74+01.32 (hereafter, G304.74). We have recently mapped 
this IRDC in the 870-//m dust continuum emission with the 
LABOCA bolometer on APEX (Paper I). The cloud was 
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found to have a integral-shaped/"hub-filament"-kind of struc- 
ture ( [Myers 2009 1, and it was resolved into twelve clumps with 
the 2D clumpf ind algorithm (Williams et al. 1994] ). Four of the 
clumps were found to be associated with MSX 8-jum point-like 
sources, indicating embedded star formation. Moreover, three 
of these M5X-bright clumps are associated with IRAS point 
sources (the IRAS sources 13037-61 12, 13039-6108, and 13042- 
6105; hereafter, IRAS13037 etc.). Eight clumps in G304.74 
were found to be dark in the MSX 8-jum image. The clump 
masses were estimated to be in the range ~ 40-220 M0 (assum- 
ing a dust temperature of 15 or 22 K), within the effective radii 
of ~ 0.3 - 0.5 pc. Unfortunately, G304.74 is not covered by the 
Spitzer infrared maps, which "only" cover the Galactic latitudes 
\b\ < r. Our 870-jum LABOCA map of G304.74 is shown in 
Fig.[T](modified from Fig. 1 of Paper I). See also Fig. 2 of Paper 
I for the MSX 8-/im image of the region overlaid with LABOCA 
contours. 

To our knowledge, besides the IRDC identification work 
by Simon et al. (2006) and our Paper I, the only studies of 
G304.74, or a clump within it, are those by Fontani et al. (2005) 
and Beltran et al. (2006). The latter authors performed SIMBA 
1.2-mm dust continuum mapping of the cloud, and identified 
eight clumps, including the M5X-bright sources (see our Paper 
I for more details). Fontani et al. (2005) carried out pointed CS 
and C'^O observations towards one source in G304.74, namely 
IRAS 13039. By employing the Galactic rotation curve of Brand 
& Blitz (1993) and the radial velocity -26.2 km s ' derived from 
CS observations, Fontani et al. (2005) determined the source 
distance of 2.44 kpc. As was recently discussed by Beuther et 
al. (2011), G304.74 constitutes a small part of the kpc-distance 
cloud complexes, which form a portion of the Coalsack dark 
cloud in the plane of the sky. 

In the present paper we will present the first molecular- 
line observations along the whole filamentary structure of 
G304.74. These data will be used to examine the basic phy- 
sical properties of the cloud, such as gas kinematics and dy- 
namical state of the clumps. The cloud's kinematic distance, 
and the distance-dependent parameters presented in Paper I, 
will also be revised. This study will also address some che- 
mical properties of the clumps. For example, we investigate 
the amount of CO depletion, which has only recently been 
studied in IRDCs ([Zhang et al. 2009] [Hernandez etal. 20111 
IMiettinen et al. 201 ll lChen et al. 20l"n i. 

The rest of this paper is structured as follows. The obser- 
vations and data-reduction procedures are described in Sect. 2. 
The direct observational results are presented in Sect. 3. In 
Sect. 4/Appendix A we describe the analysis and present 
the results of the physical and chemical properties of the 
cloud/clumps. Discussion of our results is presented in Sect. 5, 
and in Sect. 6 we summarise the main conclusions of this study. 

2. Observations and data reduction 

The observations presented in this paper were made using the 
APEX 12-m telescope at Llano de Chajnantor (Chilean Andes). 
The telescope and its performance are described in the paper by 
Glisten et al. (2006). All the submm peak positions in G304.74 
identified from our previous LABOCA 870-//m map were ob- 
served in C'^0(2 - 1), and selected positions in the southern 
part of the filament were observed in '^CO(2 - 1), SiO(5 - 4), 
and CH30H(5i; - 4^^) with the Swedish Heterodyne Facility 
Instrument (SHeFI; |Belitsky et al. 2007[ [Vassilev et al. 2008a). 
The target positions with their physical properties are listed in 
Table [1] The observed transitions and observational parameters 



are listed in Table[2| The observations took place on 18-19 and 
26-28 May 2011. As frontend for the observations, we used 
the APEX-1 receiver of the SHeFI (Vassilev et al. 2008b). The 
backend for all observations was the Fast Fourier Transfrom 
Spectrometer (FFTS: lKlein et al. 2006 J with a 1 GHz bandwidth 
divided into 8 192 channels. 

The observations were performed in the wobbler-switching 
mode with a 150" azimuthal throw (symmetric offsets) and a 
chopping rate of 0.5 Hz. The telescope pointing and focus cor- 
rections were checked by continuum scans on the planet Saturn, 
and the pointing was found to be better than ~ 5". Calibration 
was done by the chopper-wheel technique, and the output in- 
tensity scale given by the system is T^, which represents the 
antenna temperature corrected for atmospheric attenuation. The 
observed intensities were converted to the main-beam bright- 
ness temperature scale by Tmb = T'^/rjMB, where 77MB is the 
main beam efficiency. The absolute calibration uncertainty is es- 
timated to be around 10%. 

The spectra were reduced using the CLASS90 programme of 
the IRAM's GILDAS software packag^. The individual spectra 
were averaged and the resulting spectra were Hanning-smoothed 
in order to improve the signal-to-noise ratio of the data. A first- 
or third-order polynomial was applied to correct the baseline in 
the final spectra. The resulting Icr rms noise levels are ~ 40-100 
mK at the smoothed resolutions. 

The '^CO(2 - 1) and C'^0(2 - 1) rotational lines contain 
three and nine hyperfine (hf) components, respectively. We fit- 
ted these hf structures using "method hfs" of CLASS90 to derive 
the LSR velocity (vlsr) of the emission, and FWHM linewidth 
(Av). The hf-line fitting can also be used to derive the line opti- 
cal thickness, t. However, in all spectra the hf components are 
blended together, and thus the optical thickness could not be re- 
liably determined. For the rest frequencies of the '^CO(2 - 1) hf 
components we used the values from Cazzoli et al. (2004; Table 
2 therein), whereas those for C'^0(2 - 1) were taken from Ladd 
et al. (1998; Table 6 therein). 



3. Observational results 

The smoothed C'^0(2 - 1) spectra towards all target positions 
are shown in Fig. [2] The line was clearly detected towards all 
sources, but its hf components are mostly blended in all cases. 
Towards a few sources, such as SMM 2, the hf structure is just 
partially resolved, but not sufficiently well to derive the line op- 
tical thickness. Figures [3H5] show the '^CO(2 - 1), SiO(5 - 4), 
and CH30H(5t - 4^) spectra towards the observed five clumps 
in the southern part of the filament. All the '^CO(2 - 1) lines 
show a blue asymmetric profile, i.e., the blueshifted peak is 
stronger than the redshifted peak, and there is a self-absorption 
dip near the systemic velocity. The self-absorbed profiles of the 
'^CO lines indicate that the lines are optically thick. Moreover, 
the lines show non-Gaussian line-wing emission. In particu- 
lar, SMM 2 shows high-velocity blueshifted wing emission up 
to -35 km s"'. The small "absorption"-like features seen in 
the '^^CO spectra of SMM 4 and IRAS 13037 are caused by 
weak emission in the off-source reference position (off-beam). 
SiO(5 - 4) was detected only towards SMM 3, though the line is 
also very weak towards this source. We additionally detected the 
DCN(3 - 2) line towards SMM 3 in the frequency band covering 
SiO(5 - 4) (see Fig. [6]). The detected 7 = 3-2 rotational line of 
DCN is split into six hf components. In fitting this hf structure. 



' jhttp : //www ■ iram . f r/IRAMFR/GILDAS] 
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Fig. 1. LABOCA 870-yum image of the IRDC G304.74. The image is shown with linear scaling, and the colour bar indicates the flux 
density in units of Jy beam"'. The rms level is 0.03 Jy beam"' (lo"). The first contour and the separation between contours is 3cr. 
The green plus signs show the submm peak positions, which were used as the line-observation target positions in the present study. 
The white circles indicate the nominal catalogue positions of the MSX 8-yum sources. The 1-pc scale bar and the effective LABOCA 
beam HPBW (~ 20") are shown in the bottom left (modified from Paper I). 



we used the rest frequencies from the CDMS catalogue. The E- 
type methanol transition CH30H(5_i,5 -4_i_4)-E, and the A-type 
transition CH30H(5o,5 - 4o,4)-A+ at ~ 241.8 GHz were detected 
towards all the observed clumps. Moreover, towards SMM 2, a 
blend of CH30H(54,i -44,o)-A+ and CH30H(54,2 -44,i)-A- was 
detected. 

The spectral-line parameters derived from hf- and single 
Gaussian fits ai-e given in Col. (3)-(5) of Table[3] For SiO(5 -4) 
non-detections, the 3cr upper limit on the line intensity is given. 
The integrated line intensities listed in Col. (6) of Table [3] were 
computed over the velocity range given in square brackets in the 
corresponding column. The quoted uncertainties in Vlsr and Av 
are formal Icr fitting errors, whereas those in Tub and j TMsdv 
also include the 10% calibration uncertainty. In the last column 
we give the estimated peak line optical thickness, which will be 
derived in Sect. 4. 3/ Appendix A. 3. 



4. Analysis and results 

4.1. Kinematic distance of G304.74+01 .32 

The LSR velocities derived from C'^0(2 - 1) lie in the range 
Vlsr ^ [-27.4, -26.1] km s"', i.e., within 1.3 km s"'. Such ve- 
locity coherence shows that the filamentary IRDC G304.74 is a 




-40 -35 -JO ,, -,= 25 -20 -15 

Vlsr [km s ] 



Fig. 6. The smoothed DCN(3-2) spectrum towards SMM 3. The 
line is overlaid with a hf-structure fit, and the vertical lines show 
the relative velocities of individual hf components. 



coherent structure, and not just a group of physically indepen- 
dent objects seen in projection in the plane of the sky. 

The mean and median radial velocities derived from 
C'^0(2 - 1) are -26.7 and -26.5 km s ', respectively. To re- 
vise the cloud kinematic distance, we adopt the mean veloc- 
ity of -26.7 km s"', and employ the recent rotation curve of 
Reid et al. (2009). The resulting near kinematic distance, as 
expected to be appropriate for an IRDC seen in absorption, is 
d = 2.54 + 0.66 kpc. The corresponding galactocentric distance 
is about Rgc - 7.26 kpc (see Appendix A.l for further details). 
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Table 1. Source list. 



Source 


<3'2000.0 


<52000.0 


R 


M 




<n(H2)> 


MSX 8 fim 




[h:m:s] 


[°:':"] 


[pc] 


[Mq] 


[10-2 _-| 


[10* cm"''] 




SMM 1 


13 06 20.5 


-61 30 15 


0.44 ± 0.12 


121 ±92 


2.1 ± 1.2 


0.7 ± 0.5 


dark 


SMM 2 


13 06 26.9 


-61 29 39 


0.34 ± 0.09 


61 ± 46 


1.4 ± 0.8 


0.7 ± 0.5 


dark 


SMM 3 


13 06 35.8 


-61 28 53 


0.52 + 0.13 


206 ± 156 


2.3 ± 1.3 


0.7 ± 0.5 


dark 


SMM 4 


13 06 44.7 


-61 28 35 


0.47 ± 0.12 


248 ± 188 


3.0 ± 1.7 


1.1 ± 0.8 


dark 


IKAt) 13U3/-Dli/ 


13 06 49.5 


-61 Z8 U/ 


n 1 n no 
U.iZ ± U.Uo 


C/l 1 io 


1./ ± U. 1 


n o 1 n /I 
U.O ± U.4 


point 


SMM 5 


13 06 49.9 


-61 26 55 


0.32 ± 0.08 


42 ±32 


0.9 ± 0.5 


0.6 ± 0.4 


dark 


SMM 6 


13 06 51.1 


-61 27 49 


0.37 ±0.10 


109 ± 83 


2.1 ± 1.2 


1.0 ±0.8 


point 


SMM 7 


13 07 02.6 


-61 26 18 


0.37 ±0.10 


55 ±41 


1.2 ±0.7 


0.5 ± 0.4 


dark 


IRAS 13039-6108 


13 07 06.4 


-61 24 29 


0.49 ±0.13 


88 ±46 


1.2 ±0.1 


0.3 ± 0.2 


point 


SMM 8 


13 07 08.9 


-61 23 25 


0.36 ± 0.09 


67 ±50 


1.4 ±0.8 


0.7 ± 0.5 


dark 


SMM 9 


13 07 12.7 


-61 22 49 


0.39 ± 0.10 


97 ±73 


1.7 ± 1.0 


0.8 ± 0.6 


dark 


IRAS 13042-6105 


13 07 20.3 


-61 21 45 


0.33 ± 0.09 


48 ±37 


1.0 ±0.5 


0.6 ± 0.5 


point 



Notes. The clump physical properties are revised from those presented in Paper I (see Sect. 4). Columns (2) and (3) give the equatorial 
coordinates [(a, i5)2ooo.o] of the LABOCA peak position. Columns (4)-(7) list, respectively, the clump effective radius, mass, beam-averaged 
peak Ht column density, and the volume-averaged H2 number density. In the last column we indicate whether the clump appears dark or bright 
in the MSX 8-ium image. 



Table 2. Observed spectral-line transitions and observational parameters. 



Transition 


v" 


HPBW 


'7mb 


T 

± sys 


PWV 


Channel 


spacing" 




rms 




[MHz] 


["] 




[K] 


[mm] 


[kHz] 


[kms-'] 


[min] 


[mK] 


SiO(5 - 4) 


217 104.980 


28.7 


0.75 


324-360 


0.7-1.2 


122.07 


0.17 


9.5 


38^3 


DCN(3 - 2y 


217 238.61'- 


28.7 


0.75 


337 


0.8-0.9 


122.07 


0.17 


9.5 


39 


'3CO(2- 1) 


220 398.70056'' 


28.3 


0.75 


282-371 


0.8-1.0 


122.07 


0.17 


9.5 


43-46 


C"0(2- 1) 


224714.199" 


27.8 


0.75 


392-418 


0.7-1.7 


122.07 


0.16 


2.6 


81-102 


CH,OH(5_i5 -4_i4)-E 


241 767.224 


25.8 


0.75 


430-445 


0.4-0.8 


122.07 


0.15 


10.5-15.5 


39-49 



CH30H(5o.5 - 4o.4)-A+ 241791.431 

CH30H(54., -44..)-A+/-' 241 806.508 ... ... 

Notes. Columns (2)-(10) give the rest frequencies of the observed transitions (v), the APEX beamsize (HPBW) and the main beam efficiency 
('/mb) at the observed frequencies, the single-sideband system temperatures during the observations (T^ys in Tmb scale, see text), the amount of 
precipitable water vapour (PWV), channel widths (both in kHz and km s"') of the original data, the on-source integration times per position 
(fini), and the Icr rms noise at the smoothed resolution. 

The rest frequencies were taken from the Cologne Database for Molecular Spectroscopy (CDMS; 
|http : //www ■ astr o . uni -koeln . de/cdms/catalog 1 ( IMuller et al. 2005b unless otherwise stated.'*' The original channel spacings. 
The final spectra were Hanning-smoothed which divides the number of channels by two.'"' The DCN(3 - 2) transition was additionally detected 
in the band covering the SiO(5 - 4) line. The quoted frequency refers to the strongest hyperfine component F = 4 - 3.'''' The frequency was 
taken from Cazzoli et al. (2004), and it refers to the strongest hyperfine component F = 5/2 - 3/2.'"' The frequency was taken from Ladd et al. 
(1998), and it refers to the strongest hyperfine component F = 9/2- 7/2. '0 A blend of CH30H(54,i - 44,o)-A+ and CH30H(54,2 - 44,i)-A-, 
which have the same frequency. 



4.2. Revision of clump properties presented in Paper I 



We used the revised cloud distance 2.54 + 0.66 kpc to recalculate 
the distance-dependent clump parameters presented in Paper I, 
where the value d = 2.4 kpc was adopted. These include the 
clump radius (R cc d), mass (M cc d^), and volume-averaged 
H2 number density [<n(H2)) oc M/R^ ^ t/"']; see Eqs. (6)-(8) in 
Paper I. The results are presented in Cols. (4)-(7) of Table[T] The 
revised values are in the range R ~ 0.3 - 0.5 pc, M ~ 40 - 250 
Mq, A^(H2) ~ 1 - 3 X 10^2 cm-2, and (n(H2)) ~ 0.3 - 1 x lO"* 
cm For more details, see Appendix A. 2. 



Figure |7] shows the clump masses as a function of effective 
radius. A least-squares fit to the data gives log(M) - (3.08 ± 
0.23) -H (2.82 + 0.55) log(/?), with the linear correlation coefficient 
r = 0.85. For comparison, we also show the mass-radius thres- 
hold for massive star formation, M(R) = 870 M0 x (/?/pc)'^^, 
recently inferred by Kauffmann & Pillai (2010). 



4.3. Molecular column densities and fractional abundances 

Details on how the molecular column densities and fractional 
abundances were calculated are given in Appendix A. 3. In brief, 
we used Eq. ( IA.3I ) to estimate the peak optical thickness for those 
lines which were used in the column-density calculations [see 
Col. (7) of Table [3l. In all cases the lines appear to be optically 
thin(ro <K 1). Therefore, the beam-averaged column densities of 
C"0, SiO, CH3OH, and DCN were calculated using Eq. SAA[ . 
The fractional abundances of the molecules were calculated by 
dividing the molecular column density by the H2 column den- 
sity. For this purpose, the values of A^(H2) were derived from the 
LABOCA dust continuum map smoothed to the corresponding 
resolution of the line observations. The obtained column densi- 
ties and fractional abundances are listed in Table |4] 

4.4. CO depletion factors 

To estimate the amount of CO depletion in the clumps, we 
calculated the CO depletion factors, /d, following the analysis 
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-40 -35 -30 „ -^25 -20 -15 -40 -35 -30 „ -^25 -20 -15 -40 -35 -30 ,, -,,25 -20 -15 
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Fig. 2. Smoothed C'^0(2 - 1) spectra. Overlaid on the spectra are the hf-structure fits. The relative velocities of individual hf 
components are indicated with vertical lines under the spectra. 




-40 -35 -30 ,, -,,25 -20 -15 -40 -35 -30 ,, -,,25 -20 -15 

Vl5r [km , J Vl5r [km , J 



Fig. 3. Smoothed CO(2- 1) spectra overlaid with hf-structure fits. The relative velocities of individual hf components are indicated 
with vertical lines under the spectra. The red vertical line shows the systemic velocity as derived from C'^0(2 - 1). 



outlined in Appendix A.4. The derived CO depletion factors, 
0.3 - 2.3, are listed in Col. (4) of Tableg] 



4.5. Gas kinematics and internal pressure 

We used the measured C'^0(2 - 1) linewidths to calculate i) 
the non-thermal portion of the line-of-sight velocity dispersion 
(averaged over a 27'.'8 beam), and ii) the level of internal 'turbu- 
lence' (see below). The observed velocity dispersion is related 
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Fig. 4. Smoothed SiO(5 - 4) spectra. The Hne was detected only towards SMM 3 (overlaid with a Gaussian fit). Note that the 
spectrum towards SMM 3 is shown with a wider velocity range for a better illustration. 
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Fig. 5. Smoothed CH30H(5i; - 4^:) spectra overlaid with Gaussian fits. In the spectrum towards SMM 2, the CH30H(54j - 44,0)- 
and CH30H(54_2 - 44j)-A" transitions are blended. Note that the velocity range shown is different for SMM 2 for illustrative 
purposes. 



to the FWHM linewidth as CToh, = Av/ V8 In 2 ^ Av/2.355. The 
non-thermal velocity dispersion can then be calculated as (e.g., 



Myers etal. 1991 1 



(1) 



where niQUQ - 29 amu is the mass of the C'^O molecule. 
Furthermore, the level of internal turbulence is given by /mrb - 
cm/Cfi, where Cj is the one-dimensional isothermal sound speed 
(0.23 km s ' in a 15 K H2 gas with 10% He, where the mean 
molecular weight per free particle is 2.33). The values of crm 
and /tuib are given in Cols. (2) and (3) of Table |5] respectively. 
The errors in these parameters were derived by propagating the 
errors in Av and Tkii,. In Fig. |8]we plot the FWHM line widths 
as a function of clump effective radius (left panel), and the de- 
rived /turb values as a function of clump mass (right panel). The 
linewidths do not show any clear correlation with the clump size. 
Instead, they are scattered around a horizontal regression line 
(the fit is not shown), which suggests that the linewidths are more 
or less "constant" as a function of clump radius. As can be seen 



from the right panel of Fig. |8] the clumps are characterised by 
trans- to supersonic (cr^j > non-thermal motions. 

We also calculated the total (thermal +non-thermal) kinetic 
pressure within the clumps from 



^"1011 = Pt + Pm = nksT^^n + pcr^x = P(Cs + o"nt) ' (2) 

where p = 2.33mH(«(H2)) is the average gas mass density of 
the clump, and niu is the mass of a hydrogen atom. The results, 
~2-18xlO^K cm"-' (6.7 x 10'' K cm"-' on average), are Hsted 
in Col. (4) of Table|5] Note that the pressures are associated with 
considerable errors because of the assumed +5 K uncertainty in 
temperature. 



4.6. Virial analysis 

In order to examine the dynamical state of the clumps, we cal- 
culated their virial masses using the following formula, where 
the effects of external pressure, rotation, and magnetic field are 
ignored: 
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Table 3. Spectral-line parameters. 



Source 


Transition 


vlsr 


Av 




J TMBdv* 


To 








[km s"'] 


[km s"'] 


[K] 


[K km s"'] 




SMM 1 


SiO(5 - 


-4) 






< 0.13 








"C0(2 


-1) 


-26.8 ±0.01 


3.18 ±0.01 


4.92 ± 0.63(B)/4. 16 ± 0.57(R) 


15.78 ± 1.59 [-30.13,-22.65] 






C"0(2 


- 1) 


-26.3 ±0.1 


0.68 ±0.13 


0.40 ± 0.08 


0.87 ±0.14 [-28.88, -25.37] 


0.04 ± 0.02 




CH,0H(5_] 5 


- 4-i.4)-E 


-26.5 ±0.1 


1.47 ± 0.27 


0.19 ± 0.04 


0.35 ± 0.05 [-28.24, -24.57] 


0.01 ± 0.001 - 0.04 ± 0.01 




CH30H(5o.5 - 


-4o.4)-A+ 


-26. 1 ± 0. 1 


2.22 ± 0.37 


0.17 ± 0.03 


0.39 ± 0.06 [-27.57, -23.23] 


0.01 ± 0.001 - 0.03 ± 0.01 


SMM 2 


SiO(5 - 


-4) 






< 0.12 








"C0(2 


-1) 


-26.8 ± 0.03 


1.53 ±0.08 


4.57 ± 0.51(B)/2.51 ± 0.4(R) 


8.47 ±0.87 [-28.96, -23.11] 






C"0(2 


- 1) 


-27.1 ±0.1 


0.59 ±0.19 


0.53 ±0.10 


0.63 ±0.10 [-29.05, -25.66] 


0.05 ± 0.03 




CHjOHO-i.s 


- 4-i.4)-E 


-27.0 ±0.1 


0.67 ± 0.23 


0.17 ±0.03 


0.14 ±0.03 [-28.31,-26.10] 


0.01 ± 0.001 - 0.03 ± 0.01 




CH30H(5o,5 - 


^4o,4)-A+ 


-26.0 ± 0.2 


1.98 ± 0.46 


0.14 ±0.02 


0.30 ±0.06 [-28.13, -23.53] 


0.01 ±0.01 -0.03 ± 0.01 




CH,OH(54.. - 


44,)-A+'-'' 


-26.5 ± 0.04 


0.30 ± 0.80 


0.13 ±0.02 


0.17 ±0.03 [-27.02,-25.55] 




SMM 3 


Si6(5 - 


-4) 


-26.0 ± 0.5 


5.99 ± 1.06 


0.08 ± 0.05 


0.49 ± 0.09 [-32.02, -19.99] 


0.003 ±0.001 -0.07 ±0.01 




"C0(2 


-1) 


-27.7 ±0.01 


1.41 ±0.03 


8.81 ±0.93(B)/1.91 ±0.35(R) 


14.69 ± 1.48 [-29.26, -24.91] 






C"0(2 


- 1) 


-27.4 ±0.1 


0.84 ± 0.35 


0.52 ±0.11 


0.91 ±0.15 [-29.55,-25.75] 


0.05 ± 0.03 




CH30H(5_i,5 


- 4-1 4)-E 


-27.3 ±0.1 


2.70 ±0.21 


0.39 ± 0.09 


1.66±0.18[-29.75, -21.73] 


0.02 ±0.01 -0.08 ±0.01 




CH30H(5o.5 - 


-4o4)-A+ 


-26.6 ±0.2 


4.82 ± 0.60 


0.51 ±0.09 


1.85 ±0.24 [-30.75, -18.38] 


0.03 ±0.01 -0.10 ±0.01 




DCN(3 


-2) 


-27.4 ±0.1 


1.09 ±0.27 


0.13 ±0.04 


0.19 ± 0.05 [-28.76, -26.00; 92.6%] 


0.01 ±0.001 -0.11 ±0.01 


SMM 4 


SiO(5 - 


-4) 






< 0.11 








"C0(2 


-1) 


-27.6 ±0.01 


2.27 ± 0.03 


8.94 ± 0.93(B)/3.32 ± 0.42(R) 


20.04 ±2.01 [-30.13, -23.91] 






Ci'0(2 


-1) 


-27.1 ±0.03 


1.20 ±0.15 


1.33 ±0.18 


2.80 ± 0.30 [-30.55, -25.54] 


0.14 ±0.07 




CH30H(5_i,5 


- 4-,,4)-E 


-27.4 ±0.1 


1.53 ± 0.18 


0.31 ± 0.04 


0.61 ± 0.08 [-28.91, -24.40] 


0.02 ± 0.001 - 0.06 ± 0.01 




CH30H(5o,5 - 


-4o,4)-A+ 


-27.2 ±0.1 


1.47 ± 0.13 


0.40 ± 0.05 


0.60 ± 0.08 [-28.74, -25.07] 


0.02 ± 0.01 -0.08 ± 0.01 


IRAS 13037-61 12 


SiO(5 - 


-4) 






< 0.12 








"C0(2 


-1) 


-27.6 ± 0.01 


2.64 ± 0.04 


7.86 ± 0.87(B)/2.73 ± 0.46(R) 


20.66 ± 2.08 [-30.34, -23.78] 






C"0(2 


-1) 


— z /. 1 ± u.uz 


1 . lu ± U.Uo 


1 .y 1 ± u.z / 


j.yy ± [-JU.JW, -zj.zyj 


A T) a- A A 1 
U. IZ ± U.Ul 




CH30H(5_i,5 


- 4-i,4)-E 


-27.3 ± 0.03 


1.15 ±0.09 


0.53 ± 0.07 


0.73 ±0.08 [-29.41, -24.73] 


0.03 ±0.01 -0.11 ±0.01 




CH30H(5o,5 - 


-4o.4)-A+ 


-27.6 ± 0.04 


1.17±0.10 


0.43 ± 0.06 


0.70 ±0.08 [-29.41,-25.74] 


0.02 ±0.01 -0.09 ±0.01 


SMM 5 


C"0(2 


-1) 


-26.4 ± 0.04 


0.77 ±0.11 


0.81 ±0.16 


1.40±0.17 [-28.59, -25.33] 


0.08 ± 0.04 


SMM 6 


C"0(2 


-1) 


-27.0 ± 0.03 


1.66 ±0.08 


1.56 ±0.20 


3.64 ±0.38 [-30.01, -24.66] 


0.17 ±0.09 


SMM 7 


C"0(2 


-1) 


-26.3 ± 0.04 


1.36 ±0.10 


1.21 ±0.16 


2.46 ± 0.27 [-28.96, -24.45] 


0.13 ±0.07 


IRAS 13039-6108 


C"0(2 


-1) 


-26.1 ±0.04 


1.14 ±0.07 


1.12±0.14 


1.92 ±0.21 [-28.63,-24.66] 


0.07 ± 0.01 


SMM 8 


C"0(2 


-1) 


-26.2 ±0.1 


1.42 ±0.33 


0.47 ± 0.08 


1.28±0.18 [-29.17,-24.33] 


0.05 ± 0.02 


SMM 9 


C"0(2 


-1) 


-26.5 ±0.1 


2.01 ± 0.46 


0.46 ±0.10 


1.22 ±0.19 [-29.26, -25.24] 


0.05 ± 0.02 


IRAS 13042-6105 


C"0(2 


-1) 


-26.5 ±0.1 


1.00 ±0.24 


0.40 ± 0.07 


0.65 ±0.11 [-28.21,-24.66] 


0.04 ± 0.02 



Notes. Columns (3)-(8) of this table are as follows: (3) LSR velocity; (4) FWHM linewidth; (5) peak line intensity; (6) integrated intensity; (7) 
peak optical thickness (see Sect. 4.3/Appendix A.3).'''' For the '^^CO lines the intensities of the blue (B) and red (R) peak are given.**' Integrated 
intensity is derived by integrating over the velocity range indicated in square brackets. Note that the high-velocity blueshifted wing in the 
'^CO spectrum towards SMM 2 was excluded from the integrated intensity. For DCN(3 - 2), the percentage 92.6% in brackets indicates the 
contribution of hf component's intensity lying within the detected line.'''' A blend of CH30H(54_i - 44o)-A^ and CH30H(54,2 - 44j)-A". 



Table 4. Molecular column densities, fractional abundances with respect to H2, and CO depletion factors (/d). 



Source 


N{C"0) 
[10'^ cm--] 


x(C"0) 
[10-«] 


,/ 


D 


MCH3OH) 
[10" cm--] 






xCCHjOH) 
[10-'] 




SMM 1 


4.1 ± 


1.3 


3.1 4 


: 1.7 


2.2 H 


t 1.2 


2.1 : 


t 0.3^.4 + 


0.7 


1.3 


±0.8- 


-2.7 ± 


1.6 


SMM 2 


3.4 ± 


1.0 


4.1 4 


:2.3 


1.7 H 


tO.9 


1.6: 


tO.3 


-3.4 + 


0.7 


0.2 


± 0.1-0.4 ± 


0.1 


SMM 3° 


4.9 ± 


1.4 


3.0 4 


. 1.7 


2.3 H 


t 1.3 


9.94 


1.3- 


20.7 4 


:2.7 


5.8: 


t3.3- 


12.2 4 


6.9 


SMM 4 


15.1 


t3.9 


6.34 


:3.4 


1.1 H 


tO.6 


3.2: 


tO.4 


-6.7 ± 


0.9 


1.3 


±0.7- 


-2.8 ± 


1.5 


IRAS 13037-6112 


21.8 


t2.2 


21.94 


: 12.1 


0.3 H 


tO.2 


3.8: 


tO.4 


-7.9 ± 


0.9 


3.8 


±2.3- 


-7.9 ± 


4.4 


SMM 5 


7.6 4 


2.0 


13.0 


t7.2 


0.5 H 


tO.3 


















SMM 6 


19.7: 


t 5.1 


12.2 


±6.7 


0.6 H 


tO.3 


















SMM 7 


13.3: 


t3.5 


16.5 


t9.1 


0.4 H 


tO.2 


















IRAS 13039-6108 


10.5: 


t 1.1 


12.0 


t6.7 


0.6 H 


tO.3 


















SMM 8 


6.94 


1.9 


6.8 4 


3.7 


1.0 H 


tO.5 


















SMM 9 


6.6 4 


1.9 


5.94 


3.2 


1.2h 


tO.6 


















IRAS 13042-6105 


3.54 


1.0 


5.1 4 


2.8 


1.3 H 


tO.7 



















Notes. '•"> Towards SMM 3, we also determined the SiO and DCN column densities and abundances of A'(SiO) = 6.5± 1.2x 10" -2.5±0.5x 10'' 
cm-\ x(SiO) = 4.0±2.4x 10"" - 1.5 ± 0.9x 10"', andMDCN) = 6.7 ± 1.8 -48.7 ± 12.8 x 10'" cm-^ x(DCN) = 4.1 ±2.5- 30.1 ± 18.3 x 10^'". 



M„i,- = 



5 RAyj 

81n2 aG 



(3) 



where Avave is the width of the spectral line emitted by the 
molecule of mean mass fi - 2.33, and G is the gravita- 
tional constant. The parameter a = (1 - p/3)/{l - 2p/5), 
where p is the power-law index of the density profile [n(r) oc 
r'P], is a correction for deviations from uniform density 
jBertoldi & McKee 19921 hereafter, BM92; their Appendix A). 



In the clump's gravitational potential energy, the a parameter ap- 
pears as follows: 



3 GM^ 
-—a- 
5 



R 



(4) 



We note that the effect of the clump's ellipticity is neglected here 
(see BM92). We adopted the value p - 1.6 {a = 1.3) as found 
by Beuther et al. (2002) for a sample of high-mass star-forming 



7 



Miettinen: IRDC G304.74+01.32 



2.5 

2.0 

1 .5 

1 .0 

0.5 
0.2 







- 


















f 

'J — 




1 i\ 






' ¥ 













0.3 



0.4 0.5 
R [pc] 



4 T 



0.6 



0.7 



1 00 



supersonic ; 
subsonic ; 



200 300 
M \U 1 



400 



Fig. 8. Left: The FWHM linewidth of C'^0(2 - 1) versus clump effective radius. Right: The os^t/cs ratio versus clump mass. The 
dashed line represents the limit between subsonic and supersonic non-thermal motions, i.e., ctnt/cs - 1. 




1 .0 - 



-0.6 -0.5 -0.4 -0.3 -0.2 -0.1 
log(R) [pc] 

Fig. 7. Relation between mass and effective radius for the clumps 
in G304.74. The solid line shows the least-squares fit to the 
data, i.e., log(M) k 3.1 + 2.8log(R). The dotted Hne repre- 
sents the mass-radius threshold for massive star formation re- 
cendy proposed by Kauffmann & Pillai (2010), i.e., M(R) = 870 
M© X (R/pcy-^\ 



Table 5. Non-thermal velocity dispersion (ctnt), the level of in- 
ternal turbulence (/mib), and the total internal kinetic pressure 



Source 


Cr NT 


/turb 






[km 








[10^ 


SMM 1 


0.28 + 


0.06 


1.2 ± 


0.3 


2.6 


±2.0 


SMM2 


0.24 + 


0.08 


1.0 ± 


0.4 


2.2 


± 1.7 


SMM 3 


0.35 ± 


0.15 


1.5 ± 


0.7 


3.5 


±3.2 


SMM 4 


0.51 ± 


0.06 


2.2 ± 


0.5 


9.7 


±7.3 


IRAS 13037-6112 


0.46 ± 


0.03 


1.6 ± 


0.1 


6.6 


± 3.3 


SMM 5 


0.32 ± 


0.05 


1.4 ± 


0.3 


2.7 


± 1.8 


SMM 6 


0.70 ± 


0.03 


3.0 ± 


0.5 


15.3 


± 12.3 


SMM 7 


0.57 ± 


0.04 


2.5 ± 


0.5 


5.3 


±4.3 


IRAS 13039-6108 


0.47 ± 


0.03 


1.7 ± 


0.1 


2.5 


± 1.7 


SMM 8 


0.60 ± 


0.14 


2.6 ± 


0.7 


8.2 


±6.7 


SMM 9 


0.85 ± 


0.20 


3.7 ± 


1.0 


17.6 


± 15.2 


IRAS 13042-6105 


0.42 ± 


0.10 


1.8 ± 


0.5 


3.9 


±3.6 



Notes. All the parameters were calculated assuming T^^,, = 15 ± 5 K 
for all the other sources except the IRAS sources 13037 and 13039, for 
which we used the value = 22 K. 



clumps. As a function of the observed linewidth, Avobs, Av^ve is 



given by ( Fuller & Myers 1992 1 



AvL = Av^ + Av2 ^ = Av^^ + 8 In 2 X : 



mil \fi fin,o\ 



, (5) 



where Avj is the thermal linewidth, fi = 2.33, and yU^oi - 29 
for C'^O, the molecule that we have employed in the analysis. 
The virial masses are listed in Col. (2) of Table|6] The associated 
error was propagated from those of R, Avobs, and r^n. 

The virial parameters of the clumps were calculated follow- 
ing the definition by BM92, i.e.. 



2T 

W\ 



M, 



vir 



(6) 



where T is the internal kinetic energy of the clump. The value 
Q-vii- = 1 corresponds to the virial equilibrium, i.e., 2{T)+{'W) - 
0. The value o-vii < 2 corresponds to the self-gravitating limit de- 
fined by <T')-H<nV) < (e.g.. Pound & Blitz 1993 B The derived 



^ Actually, in this case, it could be more appropriate to use 
the term "equipartition" parameter instead of "virial" parameter 
dBallesteros-Paredes 20061 IHeitsch et al. 2008b . 



ffvir values are given in Col. (3) of Table |6] The uncertainty was 
derived by propagating the errors in both mass estimates. 

The derived ffyir values, ~ 0.4 - 2.8, are plotted as a function 
of clump mass in Fig.|9] Although the ofvii values are associated 
with considerable uncertainties, most of the clumps appear to 
be gravitationally bound, and some of them are near virial equi- 
librium. Moreover, there is a hint of negative correlation bet- 
ween Q-vii and M. A least-squares fit to the data points yields 
log(Q'vir) = (1.38 ± 0.63) - (0.69 ± 0.32) log(M), with the corre- 
lation coefficient r = -0.56. 



4. 7. The confining effect of external pressure 

The slope of the ayii - M relation found above, -0.69, is quite 
close to the value -2/3 ^ -0.67 characteristic of pressure- 
confined clumps (BM92). Therefore, the external pressure may 
play an important role in the dynamics of the clumps. Taking the 
kinetic energy resulting from the surface pressure on the clump, 
7"ext, into account, we can write the virial theorem in the follow- 
ing form (e.g., |McKee & Zweibel T992Pt 



^ If the magnetic energy, M, is also taken into account, the virial 
theorem is 2({T) - (Text)) + (M) + ("W) = 0. The only expansion- 
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2«r>-<rext» + <nv> = o. 



(7) 



The above energies can be written as a function of pressure as 
follows; 



and 



-3PgravVcl . 



(8) 
(9) 

(10) 



In the above formulas, Vd - ^nF? is the clump volume, Pext is 
the external pressure on the clump, and Pgrav is the gravitational 
pressure. We can write f grav as follows: 



inV _ 1/3 GM^\ 1 _ 3fl GM^ 



Substituting Eqs. dSt-lfTOl) into Eq. we get 



grav 



0. 



(11) 



(12) 



From Eq. (fTSl l we can estimate the required f ext to bring the 
clumps into virial equilibrium, which is roughly suggested by 
Fig- HI The estimated values are listed in Col. (4) of Table|6l 
We note that the values are associated with very large errors, and 
we only give the central values in the table. The clumps SMM 
\-A in the southern filament are formally characterised by ne- 
gative external pressures. Although there are large uncertainties 
in the analysis, these clumps should be either warmer or have 
some additional internal pressure acting outward, such as mag- 
netic pressure, in order to achieve virial equlibrium. The positive 
values of P^^j are found for the clumps north of SMM 4, and they 
lie in the range ~ 0.5 - 12 x 10^ K cm"^ (~ 3.3 x 10^ K cm"^ on 
average). 

4.8. Mass infall 

The detected '''CO(2 - 1) lines, most notably towards SMM 2, 
show the 'classic' signature of the infalling gas motions, i.e, 
redshifted absorption dips, or blue asymmetric profiles (e.g., 
ILucas 19761 [Myers etal.T9 96). The optically thin C'^O Unes 
peak close to the radial velocity of the self-absorption dip. 
Therefore, we can presume that the double-peak profiles are not 
caused by two distinct velocity components along the line of 
sight. Note, however, that many of the absorption features are 
slightly redshifted with respect to the systemic velocity (Fig. [3]). 
The blue asymmetry is weakest towards SMM 1, which indicates 
that in this clump there is the slowest infalling gas. We estimated 
the infall velocity, Vinf, using Eq. (9) of Myers et al. (1996): 



Vinf 



In 



)/7'dip\ 



Vied - Vblue \ 1 -I- e X T^iD/Tdip j 



(13) 



ary terms here are T and M (the angle brackets to denote time av- 
erages have been dropped for simplicity). Therefore, from the condi- 
tion Ai 2'T, we can estimate the upper limit to the magnetic field 
strength, B, when Eq. l|7j is still expected to be valid. If we neglect 
any contribution from the magnetic field outside the clump, and assume 
that the field inside the clump is uniform, we can write M = i?^ Vci/2yUo, 
where fiQ is the vacuum permeability. This leads to an upper limit of 

B « VWTVd = V6y"oXm. 



0.5 r 
0.0 r- 

^ -I.Or 
-1.5 r 
-2.0 '-_ 



_b o u n d_ 
vi ria 



1.5 2.0 
log(M) [M ; 



2.5 



Fig. 9. Relation between virial parameter and mass for the 
clumps in G304.74. The solid line is the least-squares fit to the 
data, i.e., log(avir) - 1 .4 - 0.7 log(M). The dashed line indicates 
the virial-equilibrium limit of ayir - 1, and the dash-dotted fine 
shows the limit of gravitational boundedness or ayt = 2. 

Table 6. The clump virial masses and virial parameters, and 
the estimated external pressure required to bring the clump into 
virial equilibrium. 



Source 


My 


ir 


CKvir 






[Mq] 






[10^ Kcm-3] 


SMM 1 


54 + 


20 


0.4 ± 


0.4 


-2.5" 


SMM 2 


35 ± 


16 


0.6 ± 


0.5 


-1.4" 


SMM 3 


84 ± 


54 


0.4 ± 


0.4 


-4.1" 


SMM 4 


132 d 


:44 


0.5 ± 


0.4 


-6.8" 


IRAS 13037-6112 


85 ± 


23 


1.6 ± 


0.9 


3.0 


SMM 5 


46 ± 


15 


1.1 ± 


0.9 


0.5 


SMM 6 


183 d 


:52 


1.7 ± 


1.4 


7.0 


SMM 7 


128 d 


:39 


2.3 ± 


1.9 


3.2 


IRAS 13039-6108 


137 d 


:39 


1.6 + 


0.9 


0.7 


SMM 8 


135 d 


:64 


2.0 ± 


1.8 


4.7 


SMM 9 


213 ± 


136 


2.8 ± 


2.5 


12.3 


IRAS 13042-6105 


69 ± 


32 


1.4 ± 


1.3 


1.4 



Notes. *"' The negative values suggest that stronger/additional internal 
pressure is needed to bring the clump into virial balance (see text). 



where cr is the velocity dispersion derived from appropriate op- 
tically thin line (we used the C'^O lines), Vred - Vbiue is the dif- 
ference between the radial velocities of the red and blue peaks, 
Tbd and Trd are the heights of the blue and red peaks above 
the absorption dip, and T^ip is the intensity of the absorption dip 
[see also Fig. 2 of Myers et al. (1996) for the definitions]. The 
obtained values, ~ 0.03 - 0.20 km s"', are listed in Col. (2) of 
Table |2] The uncertainty in Vinf was estimated from the error in 
the FWHM hnewidth of Ci^O(2 - 1). 

Assuming spherically symmetric infall, and that the mass 
contained in the spherical shell is dM(r) - Anr^pdr, the infall 
velocities can be used to estimate the mass infall rate as 



Mi 



inf 



dM _ dM dr 
dt ~ dr dt 



(14) 



where ri„f is the size of the infalling region (taken to be the clump 
effective radius), and p - 2.33mH(«(H2)) is again the average 
gas mass densitjQ. The derived values, ~ 2 - 36 x 10"^ M0 yr"\ 
with uncertainties estimated from the errors in Vinf , are given in 
Col. (3)ofTable|7] 



We note, however, that the density of the absorbing and infalling 
gas may be somewhat lower than the volume-average clump density. 
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Table 7. Mass-infall parameters. 



Source 


Vi„f 






Vi„f/Vff 




[kms-'] 


[IQ-^ Mq yr-'] 


[10^ yr] 




SMM 1 


0.03 ± 0.01 


2.6 ± 1.0 


45.8 ± 174 


0.02 ± 0.01 


SMM2 


0.03 ± 0.02 


1.8 ± 1.2 


33.1 ±21.4 


0.02 ± 0.02 


SMM 3 


0.17 ±0.14 


23.7 ± 19.6 


8.7 ±7.2 


0.09 ± 0.08 


SMM 4 


0.20 ± 0.05 


36.1 ±9.0 


6.9 ± 1.7 


0.09 ± 0.04 


IRAS 13037 


0.16 ±0.02 


9.8 ± 1.2 


5.5 ± 0.7 


0.13 ±0.04 



Assuming that Minf remains constant, the infall timescale can 
be calculated from 



Tinf 



M 

Mm 



(15) 



The obtained timescales, ~ 6 - 46 x 10^ yr, are listed in Col. 
of Table|7](the uncertainty is based on that in Minf). 

We also calculated the free-fall velocities of the clumps 



(4) 



Vff 



2GM 



R 



(16) 



which were then compared with the derived infall velocities. In 
the last column of Table|2]we list the Vinf/Vfr ratios. The associ- 
ated uncertainty was propagated from those of Vint , M, and R. In 
all cases, the infall velocity is much lower than the correspond- 
ing free-fall velocity. 

4.9. Fragmentation and stability analysis of the filament 

The IRDC G304.74 is highly filamentary in shape, its total mass 
is about ~ 1200 M0 (the sum of all clump masses), and it is 
resolved into twelve submm clumps at ~ 20" resolution. The 
projected linear extent of filament's long axis is about 12' or 
8.9 pc, and the filament's mass per unit length, or line mass, is 
Mine - 135 M0 pc"'. We note that this line mass may only re- 
present a lower limit, because there is likely to be a lot of mass 
on larger scales (Kainulainen et al. 201 lb). The projected sepa- 
rations between the clump submm peak positions lie in the range 
22" -115", or 0.27-1.42 peat the cloud distance (61" = 0.75 pc 
on average). The shortest separations are found near SMM 6 and 
IRAS 13037, whereas the separations are largest around SMM 7. 

To examine the fragmentation of G304.74 in more detail, we 
first calculate its thermal Jeans length from 



(17) 



where Cs was adopted to be 0.23 km s"' (Sect. 4.5), So - 
2.33mHA^(H2) is the surface density, and A^(H2) refers to the 
central column density for which we use the median value ~ 
1.4 X 10^^ cm"^ [see Col. (6) of Table[T]. The resulting thermal 
Jeans length is only ~ 0.05 pc, which is about five times less than 
the spatial resolution of our LABOCA map, and therefore much 
shorter than the observed minimum clump separation. Thermal 
Jeans instability does not appear to be responsible for the cloud 
fragmentation at the scale of clumps. 

On the other hand, we found that the filament is domi- 
nated by supersonic non- thermal motions. Therefore, it may be 
appropriate to replace c^ in Eq. ( fTTl i with the effective sound 
speed, which takes both thermal and non-thermal motions into 
account, i.e., Ces - (c^ -1- 1/3 x cr^x 3d)^^^' where (Tnt,3d 
is the three dimensional non-thermal velocity dispersion (e.g.. 



IChandrasekhar 19511 IBonazzola et al. 19871 1. Here we have as- 
sumed isotropic non-thermal motions, i.e., ctnt.sd - VSctnt, id- 
Using the median value of ctnt, id = 0.47 km s"' [see Col. (2) 
of Table |5|, the effective, or "turbulent" Jeans length becomes 
~ 0.24 pc. This is about five times larger than the thermal Jeans 
length, but still much smaller than the average clump separation 
of ~ 0.75 pc. Therefore, the observed fragmentation of the fi- 
lament cannot be attributed to non-thermal Jeans-type instability 
either. 

According to the magnetohydrodynamic "sausage"- 
instability theory, fragmentation of a self-gravitating fluid 
cylinder will lead to the formation of distinct condensations with 
almost periodic separations. This fragmentation length-scale 
corresponds to the wavelength, at which the instability grows 
the fastest (e. g.. IChandrasekhar & Fermi 1953l[Nagasawa 1987[ 
see also IJackson et al. 2 0101. If the cylinder's radius is 7?cyi, 
then, for an incompressible fluid, perturbation analysis shows 
that the above critical wavelength is /imax"'' - H^cyi- If we 
use as /?cyi the mean clump effective radius of ~ 0.4 pc, we 

obtain /tma™'' - 4.4 pc. This is clearly longer length scale than 
observed in G304.74. On the other hand, for an isothermal, 
infinitely-long gas cylinder the fastest growing mode appears at 
^max'™' - 22i/, where H ^ cj ^jAnGp^) is the isothermal scale 
height with po the central gas mass density along the axis of 
the cylinder (e.g., Nagasawa 1987). If we compute po assuming 
that the central number density is 10^ cm"^, we obtain the 
scale height H = 0.013 pc, and thus /i;^"?"™ - 0.3 pc. This is 
comparable to the efifective Jeans-length calculated above. If the 
value of H is calculated using the effective sound speed, Cetj, i.e., 
we calculate the eff'ective scale height, we get H^s - 0.03 pc, 
and /IJ^a**^™' - 0.7 pc. This wavelength is close to the observed 
average clump separation. We also note that it is comparable to 
the average width of the filament, ~ 0.8 pc, in accordance with 
the theory jNakamura et al. 19931) . 

Next, we examine whether the G304.74 filament could be 
unstable to axisymmetric perturbations or radial collapse. For 
an unmagnetised, infinite, isothermal filament the instability is 
reached if its Miine value exceeds the critical equilibrium value of 
M™^ = IcllG (e.g.. lOstriker 19641 [Inutsuka & Miyama 1992 1. 
At Tkin = 15 K, the critical value is M™^ ^ 25 M© pc"'. 
Using this value we estimate that the ratio Miine/M™^ ^ 5 for 
G304.74. Thus, G304.74 appears to be a thermally supercriti- 
cal filament susceptible to fragmentation, in agreement with the 
detected clumpy structure. If the external pressure is not neg- 
ligible, as may be the case for G304.74, then M™' is a ctually 
smaller than the above value (e.g, Fiege & Pudritz 2000a), mak- 
ing the filament even more supercritical. However, again due to 
the presence of non-thermal motions, it may be appropriate to re- 
place M™^ with the virial mass per unit length, M™^ = 2<cr^)/G, 
where (cr^) is the square of the total velocity dispersion obtained 
from Eq. © as (cr^) = Av2^,/81n2 ( |Fiege & Pudritz 2000a| ). 
Using again the value Ti^in = 15 K and the C"0 Unewidths, we 
obtain the mean velocity dispersion of (cr^)'^^ = 0.54 km s"'. 
This leads to the values M,^;;;. ^ 136 Mq pc"' and Miine/M™ ^ 
1. In this case, G304.74 as a whole would be near virial equilib- 
rium. The fragmentation timescale for a filament of radius /?cyi 
is expected to be comparable to its radial signal crossing time, 
Tcross = Rcy\/{o-y^ [scc Eq. (26) in [Fiege & Pudritz 2000b) , 



which for G304.74 is estimated to be ~ 7.2 x 10^ yr. We note 



that the value of {cr } 



2\l/2 



and thus M™ , 

line ' 



may increase if the 



cloud's gravitational potential energy is converted into gas ki- 
netic energy during the contraction of the cloud. We also note 
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that in the case of a magnetised molecular-cloud filament, the 
critical mass per unit length is unlikely to differ from that of 
an unmagnetised filament by more than a factor of order unity 
(|Fiege & Pudritz 2000a|b). 



5. Discussion 

5.1. Fragmentation and dynamical state of G304.74+01 .32 

The average projected clump separation in G304.74 is about 0.75 
pc. As shown above, this value can be explained if the frag- 
mentation of the parent filament is caused by "sausage"-type 
(or pinch) instability, where the dominant unstable mode ap- 
pears at the wavelength depending on the effective sound speed. 
Such instability could have been triggered by, e.g., increased ra- 
dial inward force or radial contraction leading to the increase 
of the tension force of the azimuthal magnetic field. This con- 
forms to the results of Jackson et al. (2010), who speculated that 
the "sausage"-type instability is the main physical driver of frag- 
mentation of filamentary IRDCs. Henning et al. (2010) and Kang 
et al. (2011) recently found that the mean clump separation is 0.9 
pc in both the filamentary IRDC GO 1 1 . 1 1 -0. 1 2 or the Snake, and 
IRDC G049.40-00.01, respectively. This resembles the situation 
in G304.74, and could therefore suggest a similar fragmenta- 
tion process. In Paper I we discussed the possibility that the ori- 
gin of the G304.74 filament, and its fragmentation into clumps 
is caused by supersonic turbulent flows. Even the formation of 
the cloud could have been caused by converging turbulent flows, 
the role of random turbulence in fragmenting the filament ap- 
pears unclear in the light of the present results. Nevertheless, the 
new results support the importance of non-thermal motions in 
G304.74. 

The slope of the a^„-M relation we found is roughly consis- 
tent with the prediction of BM92 for clumps confined by ambi- 
ent pressure from their surrounding medium, i.e., ffvir °^ M^^^^ . 
This slope can be understood as follows. On one hand, the clump 
mass was found to be roughly proportional to (Fig. |7]l. On 
the other hand, the linewidth-radius 'relation' shown in Fig. [8] 
is roughly flat, i.e., Av ^constant. Therefore, from Eqs. (|3]l and 
(|6]i we see that ofvii- oc RAv^/M cc M^I^M'^ cc M'^^^. Although 
our clumps, or most of them, appear to be gravitationally bound 
(ttvir ^ 2), it is possible that the ambient pressure (still) plays 
a non-negligible role in the clump dynamics. Note that for real 
pressure confined clumps, o-vii s> 1, and they are not gravitation- 
ally bound (BM92). 

Could the external pressure be caused by turbulent flow mo- 
tions outside the cloud ? To estimate the turbulent ram pres- 
sure, Pram - pcT^x' ^c usc the average '^CO(2 - 1) linewidth 
of ~ 2.2 km s ' derived for five of our clumps, and assume 
the density of the '^CO emitting gas to be lO^* cm"^ (fol- 
lowing ' Lada et al. 20081 ). If T^a 

is in the range 10—20 K, the 
corresponding non-thermal velocity dispersion is ctnt - 0.93 
km s"', which yields the value Prum/ks ~ 3 x 10^ K cm"-'. 
This is only about two times lower than the average internal ki- 
netic pressure within the clumps, and very close to the average 
external pressure required for virial equilibrium (Pram - Pext' 
Sect. 4.7). Therefore, this supports the scenario of the impor- 
tance of turbulence presented in Paper I. The obtained Pram/^B 
value is also consistent with intercloud pressures of 10"* - lO*" 
K cm"^ deduced from theoretical estimates and different obser- 
vations (see Field et al. 201 1 and references therein). However, 
we note that the large-scale turbulent motions are anisotropic, 
and therefore the associated pressure may not be able to confine 
the cloud, but could rather distort and/or even disrupt it (e.g.. 



IBallesteros-Paredes 2006] l. The finding that Mn^^/M™^ ^ 1 for 
the G304.74 filament suggests that it is close to virial equili- 
brium. The virial balance may be difficult to explain in the 
context of turbulent fragmentation discussed in Paper I, be- 
cause turbulent motions are expected to induce a flux of mass, 
momentum, and energy between the cloud and its surround- 
ing medium (I Ballesteros-Paredes 2006). Also, the detection of 
large-scale infall motions towards G304.74 suggests it to be out- 
of-equilibrium. For comparison, the recent Herschel study of 
the Aquila rift and Polaris Flare regions by Andre et al. (2010) 
showed that practically all filaments with Miine/M™^ > 1 were 
associated with prestellar cores and/or embedded protostars. On 
the other hand, filaments with Miine/Mj 



line 



1 were generally 

found to lack such star-formation signatures. More observations 
with higher resolution would be needed to better characterise the 
star-formation activity in G304.74, and to compare it with fila- 
mentary structures revealed by Herschel. 

If the Miine/M™^ ratio is indeed ^ 1, then from Eq. (11) of 
Fiege & Pudritz (2000a), and the average ratio of the external 
to internal pressure of (Pext)/(^kin) - 0.49, we estimate that the 
ratio between the total magnetic and kinetic energies per unit 
length is MIYW\ - 0.49. This value is > 0, and therefore sug- 
gests that i) the overall effect of the magnetic field on the cloud 
is to provide support, and ii) the magnetic field is poloidally do- 
minated. This is in contrast with the finding of Fiege & Pudritz 
(2000a), i.e., that the magnetic fields in the filamentary clouds 
they analysed are likely to be helical and toroidally dominated. 
However, Fiege et al. (2004) analysed the Snake IRDC, which 
is more similar to G304.74, and concluded that it is likely do- 
minated by the poloidal component of the magnetic field (or 
is magnetically neutral), and is therefore an excellent candi- 
date for a magnetically supported filament. This may also be 
the case for the G304.74 filament. Interestingly, Fiege & Pudritz 
(2000b) found that models with purely toroidal fields are sta- 
ble against "sausage" modes, whereas models with a poloidal 
field component may be unstable. We note, however, that the 
observed fragmentation of the filament strongly suggests that 
the magnetic field is not purely poloidal (or toroidal), because 
otherwise the fragmentation would have presumably been sta- 



bilised (Fiege & Pudritz 2000b). If the filament fragmentation is 



caused by sausage-type instability, some toroidal magnetic field 
is probably needed, because magnetic pressure associated with 
the poloidal field would resist the "squeezing" of the cloud. 

On the other hand, the ambient pressure can also arise from 
the cloud's self-gravity. This was found to be the case in the Pipe 
Nebula by Lada et al. (2008), where all the cores within the cloud 
where found to have avii > 1, and avh <^ M^^-^^. The pressure 
due to the weight of the cloud is 



37r. 
20 



(18) 



where S - M/nR^ is the mean mass surface density of the cloud, 
and <pG is a dimensionless correction factor to account for the 
non-spherical geometry of a cloud [BM92; cf. our Eq. (fTTl ll. 
Following BM92, we use as the mean effective cloud radius 
the value R ^ -y/0^8^^"8^9pc2 ^ 2.7 pc, where 0.8 pc is the 
average width of the filament. The corresponding surface den- 
sity is E ^ 1 g cm"^. The aspect ratio of the filament is about 
y - 8.9pc/0.8pc ~ 11 (prolate), and therefore we adopt the 
value - 2 (see Appendix B of BM92). Thus, we obtain the 
value Pel ~ 5.5 X 10'* K cm . This is about an order of mag- 
nitude lower than the average internal kinetic pressure within 
the clumps or the average P^^, value, suggesting that the cloud's 
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weight is not contributing significantly to the overall confining 
pressure. We note, however, that Kainulainen et al. (201 lb) re- 
cently found that the masses of IRDCs traced by the thermal dust 
emission are generally only < 10 - 20% of the masses traced by 
near-infrared and CO data. Therefore, the pressure due to the 
weight of the cloud itself, which is mostly caused by the mate- 
rial at low column densities, is likely to be higher than estimated 
above. 

Some of the clumps in the southern end of the filament were 
found to be characterised by negative external pressures in the 
case of virial equilibrium (or 'energy equipartition' ; Sect. 4.7), 
suggesting that they should be either warmer or have some ad- 
ditional internal pressure. Interestingly, this conforms to the re- 
sults obtained in Paper 1, where we used the dust opacities at 
8 fim and 870 jum to estimate the dust temperatures towards the 
MlR-dark submm peaks. The clumps SMM 1 and 2 in the south- 
ern tip of the filament were estimated to have temperatures of 
~ 30 K, but the values were associated with very large relative 
errors of ~ 60 - 70% (therefore we have adopted Tdust = 15 
K). The masses (mean densities) of SMM 1 and 2 would be 
~ 44 Mq (0.2 X 10"* cm-3) and 22 M© (0.3 x 10^ cm"^), re- 
spectively, if Tdust = 30 K. The non-thermal velocity dispersions 
would be 0.27 and 0.23 km s"', and the total kinetic pressures 
about 1.2 X 10^ K cm"^ and 1.6 x 10^ K cm"^, respectively. 
Consequently, the external pressures required to virial balance 
would become positive, i.e., fg^(/A;B ~5.2xl0^and~ 1.1x10^ 
K cm""*, for SMM 1 and 2, respectively. These are in reasonable 
agreement with the values estimated towards other clumps 
[see Col. (4) of Table|6J. The present results seem to support the 
possibility of north-to-south temperature gradient of ~ 15 — > 30 
K in G304.74. 



5.2. Clump kinematics, infall motions, and prospects for 
massive star formation 

The clumps in G304.74 appear to be characterised by trans- 
to supersonic internal non-thermal motions. According to 
the Turbulent Core Model by McKee & Tan (2003), the 
clump (or core) has to be supersonically turbulent in order 
to achieve a high mass accretion rate and be able to form 
massive stars. The presence of such turbulence is also ex- 
pected to lead to the formation of substructure within the 
clump. Higher-resolution observations would be needed to 
find out whether the clumps are fragmented into smaller 
cores. For example, interferometric observations of a sample 
of clumps associated with IRDCs have revealed the presence 
of embedded cores within them (e.g., IRathborne et al. 20081 
IHennemann et al. 2 009; Wa ng et al. ZoTT) . 

The clump SMM 6, which has the second largest /mrb value 
(~ 3), is associated with two MSX point sources, suggesting that 
it indeed hosts subfragments, and that they have already been 
collapsed. On the other hand, as shown in Fig. |7] none of the 
clumps appear to clearly fullfil the mass-radius threshold for 
massive star formation proposed by KaufFmann & Pillai (2010). 
Only within the uncertainties, the two most massive clumps, 
SMM 3 and 4, could exceed this Umit. 

The supersonic "turbulent" motions in the clumps 
may actually be due to gravitational infall motions (see 
IHeitsch et al. 20091 and references therein). Evidence of large- 
scale infall motions in '^CO lines was indeed found towards the 
southern clumps (see further discussion below). If the whole 
filament is in a process of hierarchical gravitational collapse 
with clumps collapsing locally, a complex velocity pattern is 



expected to be created. Such "turbulence" is not able to resist 
gravitational collapse; instead, it is generated from the collapse 
itself (Ballesteros-Paredes et al. 201 1 ). 

Radiative transfer models have shown that if the clump 
centre is collapsing faster than the outer layers, the observed 
spectral line should have two peaks with the blueshifted 
peak being stronger than the redshifted one (Zhou et al. 19931 
Myers et al. 1996) . The blue peak originates from the rear side 
of the clump, whereas the red peak comes from the front part 
of the clump. The central dip between the two peaks is caused 
by the envelope layer (see also lEvans 19991 1. Such line profiles 
were seen towards SMM 1-4 and IRAS 13037, and therefore 
they are likely undergoing large-scale inward motions. We de- 
rived the infall velocites and mass infall rates of ~ 0.03 - 0.20 
km s"' and ~ 2 - 36 x 10"^' Mq yr"\ respectively. These are 
comparable to the values 0.14 km s"' and 3 x 10"^ Mq yr"' 
found by Hennemann et al. (2009) towards J 18364 SMM 1 
South, which is a massive clump within an IRDC. Chen et al. 
(2010) measured higher infall velocities and rates of ~ 0.7 - 
4.6 km s"' and ~ 1 - 28 x lO^^* Mq yr"', respectively, for 
their sample of MYSOs associated with extended green objects 
(EGOs) and IRDCs. Typically, the m ass infaU rates in MYSOs 
are ~ 10 ^^ - 10;^ Mq yr ' (e.g., IKlaassen & Wilson 20071 
I Zapata et al. 2008[ IBeltran et al. 201 11 1. We note that the infall 
timescales we derived, ~ 6 - 46 x 10^ yr, are quite long com- 
pared to the characteristic timescale of ~ lO'' yr for high-mass 
star formation (e.g.. IDavies et al. 201 II and references therein). 
The infall rate, however, may not stay constant during the star- 
formation process, which could also explain some of the large 
observed variations in values quoted above. 

We found that the infall velocities are one to two or- 
ders of magnitude lower than what would be expected from 
a free-fall collapse. This suggests that the large-scale col- 
lapse is not gravitationally dominated, but it could be retarded 
by, e.g., magnetic fields. The predominantly poloidal mag- 
netic field could resist motions perpendicular to the poloidal 
field lines, such as appears to be the case in the Snake IRDC 
( |Fiege & Pudritz 2000bl [Fiege et al. 2004). We also emphasise 
that the derived infaU parameters refer to the infall of the sur- 
rounding gas onto a (pre-)protocluster, not onto individual stars. 
For example, the timescale for (disk) accretion from the infalling 
envelope may be substantially shorter than the infall timescale. 
The current data do not allow us to decide whether the ob- 
served large-scale mass infall could feed competitive accretion 
in the clumps/protoclusters dBonnell & Bate 20061 1. or be re- 
lated to monolithic collapse leading to massive star formation 
(IMcKee & Tan 2003i l. On one hand, the formation of high-mass 
stars in G304.74 seems uncertain on the basis of the mass-radius 
threshold proposed by Kauffmann & Pillai (2010). On the other 
hand, the high luminosites of the IRAS sources 13037 and 13039, 
~ 1.5 - 2 X lO-' Lq (Paper I), shows that at least intermediate- 
mass star- and/or cluster formation is taking place in this IRDC. 



5.3. The observed cliemical properties of the clumps 

5.3.1 . CO depletion and the SiO detection towards SIVIM 3 

The CO depletion factors we derived are in the range fu ~ 
0.3 - 2.3. Therefore, CO molecules do not appear to be signi- 
ficantly depleted at the scale of clumps in G304.74, which con- 
forms to the low average H2 densities. Higher-resolution obser- 
vations would be needed to investigate whether the clumps con- 
tain smaller and more CO-depleted cores. The values /d < 1 
could be due to too small a canonical CO abundance used, or 
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due to contamination of non-depleted gas along the line of sight 
(see also IMiettinen et al. 201 H and references therein). The lat- 
ter possibility, however, may be less important because the cloud 
is located above the galactic plane, and thus relatively free of 
contaminating emission along the line of sight. Of course, uncer- 
tainties in A^(H2), particularly due to the dust opacity, propagate 
into jc(C'^O) and /o. 

It is not surprising that we found no evidence of CO freeze- 
out towai-ds the IRAS sources and SMM 6, which ai-e all MIR 
bright and very likely associated with ongoing star formation. 
The physical processes associated with star-formation activity, 
such as outflows, heating, etc., are expected to release CO from 
the icy grain mantles back into the gas phase. Consistently, the 
largest fu values are seen towards the MIR-dark clumps SMM 
1 and 3. However, the detection of SiO emission in SMM 3 
suggests the presence of outflows, because gas-phase SiO is ex- 
pected to be produced when shocks disrupt the dust grains (e.g., 
Hartquist et al. 19801 jSchilke et al. 1997] ). This is in agreement 



with the finding of Sakai et al. (2010; henceforth, SSHIO), i.e., 
that the fraction of shocked gas (as traced by SiO) is higher in the 
MSX-dark sources than in the MS'X-bright sources. This could 
indicate that the MIR-dark clumps are associated with newly 
formed shocks giving rise to stronger SiO emission, whereas the 
weaker emission from the MIR-bright clumps would be due to 
shocks being older, and thus milder. We note that the highest es- 
timated SiO column density value in SMM 3, ~ 2.5 x 10'^ cm"^, 
is comparable to the values ~ 1.1 - 2.6 x lO'"* cm"^ found by 
SSHIO towards their MSX-dark clumps. 

For comparison, Hernandez et al. (2011) found depletion 
factors up to /d ~ 5 in the highly filamentary IRDC G035.39- 
00.33. Miettinen et al. (201 1) derived the values /d ~ 0.6 - 2.7 
towards a sample of clumps within different IRDCs, which are 
very similar to those obtained in the present work with the same 
telescope/angular resolution. Higher values of /d ~ 6 - 19 were 
recently found by Chen et al. (2011) for their sample of high- 
mass star-forming clumps, including sources within IRDCs, u- 
sing the 10-m SMT observations of C'^O at 34" resolution. This 
shows that strong CO depletion is possible even at the scale of 
clumps. We note that Hernandez et al. (201 1) used the CO deple- 
tion timescale to constrain the lifetime of G035. 39-00. 33 (or the 
part of the cloud showing the strongest depletion). Such inves- 
tigations are useful considering the cloud formation timescales 
and mechanisms. Unfortunately, the fact that we found no clear 
evidence of CO depletion in G304.74 does not allow us to per- 
form similar analysis. 

5.3.2. CH3OH 

The fractional CH3OH abundances in the observed clumps were 
estimated to lie in the range ~ 0.2 - 12 x 10"^. The lowest values 
are found towards the MIR-dark clump SMM 2, and the highest 
ones towards SMM 3 and IRAS 13037. 

The pure gas-phase chemical models suggest that CH3OH 
is formed from methylium cation through the radiative associ- 
ation reaction CHj + H2O CH3OH2 + hv, followed by the 
electron recombination reaction CH3OH2 H- e — > CH3OH + H. 
The latter reaction is, however, too inefficient to explain the ob- 
served CH3OH abundances in star-forming regions alone (see, 
e.g.. IWirstrom et al. 201 H and references therein). For example, 
the CH3OH observations by van der Tak et al. (2000) towards 
a sample of 13 massive young stellar objects (MYSOs) sug- 
gested three types of abundance profiles: the coldest sources 
have CH3OH abundances of ~ 10 warmer sources are cha- 
racterised by a range of abundances from ~ 10"^ to ~ 10"^, 



whereas the hot cores have a:(CH30H) ~ 10"^. These high abun- 
dances indicate that solid CH3OH must have been evaporated 
from the icy grain mantles into the gas phase. In the grain- 
surface chemistry scheme, CH3OH is formed through succes- 
sive hydrogenation of CO on dust grains (CO — > HCO — > 
H2CO — > CH2OH — > CH3OH), and subsequent grain man- 
tle evaporation is then predicted to cause the gas-phase abun- 
dance of CH3OH to be > lO '^ (e.g., Ivan der Tak et al. 20001 
Charnley et al. 2004[ IGarrod et al. 20081) . Besides grain heat- 
ing, methanol can be released from grain mantles by moderate 
shocks (e.g., IBachiller & Perez Gutierrez 19971 and references 
therein). 

What comes to the observed methanol abundances in the 
present work, the lowest values, which are found in SMM 2, do 
not appear to necessarily require mantle evaporation (see below). 
The detection of high-velocity '^CO wing emission in SMM 2, 
however, suggests the presence of outflow activity. The rest of 
the x(CH30H) values, particularly in SMM 3 and IRAS 13037, 
appear to require that CH3OH has been injected from dust grains 
into the gas. In SMM 3 the high CH3OH abundance together 
with the SiO detection and '^CO wing emission, are likely re- 
lated to outflow activity. The latter wing-emission signature was 
also found towards IRAS 13037, but the SED of this source 
suggests the presence of hot dust (Paper I). Therefore, in ad- 
dition to grain processing in shocks, the central YSO(s) may 
have raised the dust temperature above 100 K in IRAS 13037, 
which is the approximate sublimation temperature of CH3OH 
(e.g., IGarrod & Herbst 20061 ). Interestingly, besides the highest 
CH3OH abundance, the strongest CO depletion was found to- 
wards SMM 3. This could be a sign of enhanced methanol pro- 
duction in SMM 3 due to CO freeze-out (e.g., Whittet et al. 2011] 
and references therein); we still see a remnant of stronger CO de- 
pletion when CO hydrogenation on dust grains was arised, but 
also the product of this process, i.e., the gas-phase CH3OH. 

The CH3OH observations towards IRDCs have also been 
done by other authors. Beuther & Sridharan (2007) studied a 
sample of massive clumps within IRDCs, and found CH3OH 
column densities and abundances of 1.7 - 37.2 x lO'"* cm"^ and 
0.4 - 10.7 X 10"'", respectively. Although their column densi- 
ties are very similar to our values, our abundances are generally 
higher As discussed by Beuther & Sridharan (2007), their a- 
verage CH3OH abundance of 4.3x 10 '° is similar to the CH3OH 
abundances of 6 x 10"'" found in the low-mass starless cores 
L1498 and L1517B by Tafalla et al. (2006). There are no inter- 
nal heating/outflows in starless cores, so grain processing is not 
expected in such sources. Sakai et al. (2010) derived CH3OH 
column densities in the range < 2.9 - 24.8 x 10'^ cm"^ towards 
a sample of 20 massive clumps within IRDCs, which are much 
higher than those found in the present work. Sakai et al. (2010) 
pointed out that their much higher A^(CH30H) values compared 
to those of Beuther & Sridharan (2007) could be due to the 
source selection criteria; SSHIO studied MSX-Aaxk clumps with 
signatures of star-formation activity. Based on the observed cor- 
relation between the CH3OH and SiO abundances, SSHIO con- 
cluded that the origin of CH3OH in the MS'X-dark clumps is re- 
lated to shocks. Gomez et al. (201 1) found CH3OH abundances 
of~4x 10'' -3 X 10"^ in the IRDC core Gl 1.1 1-0.12 PI, and 
also suggested that outflow desorption is responsible in releas- 
ing CH3OH from ice mantles. With comparable CH3OH abun- 
dances, this could also be the case in the MIR-dark clump SMM 
3. 
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5.3.3. DCN 

We detected the 7 = 3 - 2 transition of DCN towards SMM 3. 
To our knowledge, this is the first reported detection of DCN 
in an IRDC [Sakai et al. (2012) recently detected the DNC 
isotopomer towards different IRDC clumps]. The first detec- 
tion of DCN in interstellar medium was made by Jefferts et 
al. (1973) towai-ds the Orion Nebula. Mangum et al. (1991) 
mapped the Orion-KL region in DCN, and they suggested that 
most of the gas-phase DCN could have been released from icy 
grain mantles. Both the DCN and DNC molecules are the ma- 
jor products of the association reaction D -i- CN on grain sur- 
faces (IHiraoka et al. 20061) . and DCN can be thermally desorbed 
at ~ 100 K (e.g.. Albertsson et al. 201 1). In SMM 3 we may see 
shock-originated DCN emission, which conforms to the results 
discussed above. 



6. Summary and conclusions 

We have carried out a molecular-line study of the IRDC G304.74 
using the APEX telescope. All the clumps along the filamentary 
cloud were observed in C'^0(2- 1), and selected positions in the 
souther part of the filament were observed in the rotational lines 
of SiO, '^^CO, and CH3OH. The main results are summarised as 
follows: 

1. The ~ 9 pc long filamentary IRDC G304.74 is spatially co- 
herent structure as indicated by the uniform C'^0(2 - 1) ra- 
dial velocities along the cloud. 

2. The fragmentation of the filament into clumps appears to 
be caused by "sausage"-type fluid instability, in agreement 
with results from other IRDCs, such as the "Nessie" Nebula 
dJackson et al. 2010 ). At smaller scales, however, clumps 
may have fragmented into smaller cores because of, e.g., 
Jeans-type instability. Studying this possibility would re- 
quire higher-resolution observations. 

3. Most of the clumps appear to be gravitationally bound 
(avii ^ 2), although we found evidence that the external pres- 
sure may play an important role in their overall dynamics. 
The ambient pressure could be due to the turbulent ram pres- 
sure, in agreement with the hypothesis that the formation of 
the whole filament is caused by converging supersonic tur- 
bulent flows (Paper 1). 

4. The whole filament might be close to virial balance, as sug- 
gested by the analysis of the virial mass per unit length. 
Radial collapse could be retarded by the poloidal magnetic 
field component. Clump-by-clump virial analysis suggests 
that some of the clumps in the southern part of the filament 
may be warmer than the rest of the clumps. This is consis- 
tent with our analysis of dust opacity ratios at 8 and 870 jjm 
in Paper I, which suggested dust temperatures of ~ 30 K in 
SMM 1 and 2, although with considerable uncertainties. 

5. The clumps show trans- to supersonic non-thermal internal 
motions. Moreover, the observed clumps in the southern part 
of the filament show large-scale infall motion signatures in 
the '^CO(2-l) line. The infall velocities and mass infall rates 
were estimated to be ~ 0.03 -0.20 km s"' and - 2-36x10"^ 
Mq yr"', respectively. 

6. None of the clumps clearly fulfill the mass-radius limit of 
massive star formation recently proposed by Kauffmann & 
Pillai (2010). The clumps may "only" form stellar clusters 
and/or intermediate-mass stars. Two IRAS sources in the 
cloud have luminosities ~ 1.5 - 2 x 10^ Lq (Paper I), which 
supports the latter scenario. 



7. The CO molecules do not appear to be significantly depleted 
in the clumps (/d < 2). The star-formation activity, such as 
outflows, may have released CO from the icy grain mantles 
into the gas phase. Some of the '^CO(2 - 1) lines also show 
broad wings of emission, suggesting the presence of outflow- 
ing gas. 

8. We estimate that the fractional CH3OH abundances in the 
clumps are around 2 x 10"'" - 1 x 10"**. No evidence of hot- 
core type methanol abundance of ~ 10"^ was found among 
the observed five clumps. Outflow activity appears to be a 
likely explanation for the observed CH3OH abundance in the 
MlR-dark clump SMM 3. The presence of outflows in SMM 
3 is also indicated by the detection of SiO and DCN emission 
towards this source. 
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Appendix A: Derivation of the physical and 
chemical properties of the cloud/clumps 

This appendix presents details on the analysis and formulas dis- 
cussed in Sects. 4. 1 - 4.4. 



A.1. Kinematic distance 

To calculate the cloud kinematic distance, we adopted the mean 
C'^0(2 - 1) radial velocity of -26.7 km s We also employed 
the recent rotation curve of Reid et al. (2009), which is based on 
measurements of trigonometric parallaxes and proper motions 
of masers in high-mass star-forming regions. The best-fit rota- 
tion parameters of Reid et al. (2009) are (0o, /?o)=(254 km s"', 
8.4 kpc), where Rq is the solar galactocentric distance, and 0o 
is the rotation velocity at /?(>. The resulting near kinematic dis- 
tance, as expected to be appropriate for an IRDC seen in absorp- 
tion, is c/ = 2.54 ± 0.66 kpc. The corresponding galactocentric 
distance is about Rcc - 7.26 kpc. Actually, besides being an 
IRDC, the galactic latitude of G304.74, b = -h1°32, supports the 
idea that it lies at the near distance, because then the cloud is 
~ 59 pc above the galactic plane. In the case of the far kine- 
matic distance, which is ~ 7.05 kpc, it would be suspiciously 
high above the galactic plane, i.e., ~ 162 pc. This would be more 
than twice the scale height of the molecular disk, zi/2 ~ 70 pc 
(IBronfman et al. 19881 see also lFontani et al. 2005l l. 



A.2. Revision of clump properties presented in Paper I 

The revised cloud distance was used to recalculate the distance- 
dependent clump parameters presented in Paper I, i.e., the 
clump radius (R cc d), mass (M cc cfi), and volume-averaged 
H2 number density [<n(H2)> oc M/R^ oc c/^']; see Eqs. (6)-(8) 
in Paper I. Previously, we assumed a uniform dust temperature 
of Tdust - 15 K for all the other sources except IRAS 13037 
and IRAS 13039, for which we used the value ^ 22 K derived 
from their spectral energy distributions (SEDs; see Fig. 5 of 
Paper I). Here, we adopt these same temperature values but 
the assumed 15 K temperature is now assumed to have a +5 K 
uncertainty. This is expected to be a reasonable choice, because 
previous molecular-line observations of clumps within IRDCs 
have shown the typical gas kinetic temperature to he in the 
ran ge Tkin ~ 10 - 20 K ( |Carey et al. 1998} [Teyssier et al. 2002 
Srid haran et al. 20051 



Pillai et al. 20"06r ISakai et al. 2008 



IDevine etal. 20111 



Ragan et al. 20TT] i. 

where 



cm 



Zhang et al. 2011 

Moreover, at high densities of n(H2) > 3 x 10"* 
collisional coupling between the gas and dust becomes efficient, 
the gas and dust temperatures are expected to be similar (e.g., 
IGalli et al. 2002l l. For comparison, Nguyen Luong et al. (2011) 
recently derived a mean dust temperature of ~ 14.5 K for the 
clumps/cores in the filamentai-y IRDC G035.39-00.33. To esti- 
mate the errors in R, M, and (m(H2)), we took the temperature 
uncertainty, and the distance error of +0.66 kpc into account. 
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The ±5-K temperature error was also used to estimate the 
uncertainty in the H2 column density. 

A.3. Molecular column densities and fractional abundances 

The spectral line's intensity can be expressed as the main-beam 
brightness temperature, which is given by 

Tub = ^ \f(T,,) - FinJ (1 - e'^^) , (A.l) 

where h is the Planck constant, is the Boltzmann constant, Tex 
is the line excitation temperature, Jbg is the background tempe- 
rature (taken to be the cosmic background radiation temperature 
of 2.725 K), Ty is the line optical thickness, and the function 
F{T) is identical to the (average) photon occupation number: 



F(T) = 



1 



,hv/ksT _ I 



(A.2) 



We note that the often used Rayleigh- Jeans equivalent tempera- 
ture is given by J{T) = hv/kB x F(T). From Eq. dA.lb we can 
solve for the peak optical thickness: 



To = -ln<^l- 



' MB 



j;[f(T,,)-F(T^^ 



(A.3) 



Equation iA.3\ was used to estimate the value of tq for those 
lines which were used in the column-density calculations [see 
Col. (7) of Table[3|. The assumptions made concerning the value 
of Tex are described below. 

The beam-averaged column densities, A^(mol), of C'^O, SiO, 
CH3OH, and DCN were calculated under the assumption of lo- 
cal thermodynamic equilibrium (LTE), and assuming optically 
thin emission. The latter assumption seems to be valid for our 
lines. In this case, A^(mol) is directly proportional to the inte- 
grated line intensity as 



A'(mol) 



3^Beo 1 Z,ot(rex) 

f(rbg) 



27r2 v/jIs 



gKgl 1 - 



rMfidv, (A.4) 



where eo is the vacuum permittivity, //ei is the permanent electric 
dipole moment, S is the line strength, Zrot is the rotational par- 
tition function, gx is the ^T-level degeneracy, gi is the reduced 
nuclear spin degeneracy (see, e.g., [Turner 1991 i and is the 
energy of the upper- transition state. We note that in the deriva- 
tion of Eq. ( IA.4I ). the electric dipole moment matrix element is 
defined to be - l^liS /gu, where gu - 2J + I is the rotational 
degeneracy of the upper state dTownes & Schawlow 1975] ). Note 
that for linear molecules, ^a: = ^/ = 1 for all levels, whereas 
for A-type CH3OH, gK = I and gi = 2, and for Zs-type 
CH3OH, gK = 2 and gj = 1 ( ITurner 199li We derived 
the CH3OH column densities from the integrated intensity of 
both the CH30H(5o,5 - 4o,4)-A+ and CH30H(5-i,5 - 4_i,4)-E 
lines, and the results were found to be comparable to each 
other. In this paper we only give the A^(CH30H) values de- 
rived from the above A^-transition. Note that we have detected 
too few lines to employ the population diagram analysis to de- 
rive the rotational temperature and column density of CH3OH 
([Goldsmith & Langer 1999 ). 

For the linear molecules C'^O, SiO, and DCN, we approxi- 
mated the partition function as 



kB Tex 1 
Zrot(^ex) — TT; T ' 

hB 3 



(A.5) 



where B is the rotational constant. The above expression is ap- 
propriate for heteropolar molecules at the high temperature limit 
(/iB/feeT^ex ^ !)■ The partition function of methanol, which is 
a slightly asymmetric rotor with one internal rotor and A- and 
/i-symmetry states, was calculated as 



Z,-ot(rex) = Z,o,(A)+Zro,(£) = 2 



7r{kBT,,y 
h^ABC 



(A.6) 



where A, B, and C are the rotational constants 
(ITurner 19911 ). The spectroscopic properties of the observed 
molecules/transitions are given in Table I A. ll 

The unknown parameter in the above optical-thickness and 
column-density calculations is the excitation temperature, Tex. 
In the case of C'^O we assumed that the line is thermalised at 
the clump temperature, and adopted the value Tex - 15 + 5 K 
for all the clumps except the IRAS sources 13037 and 13039, for 
which the value Tgx = 22 K was used (Appendix A.2). This is 
likely to be true for our clumps, because even their average H2 
densities are comparable to the critical density of C'^0(2 - 1) 
[9.5 X 10^ cm-^ see Col. (2) of Table|AD. Fontani et al. (2005) 
actually derived the C'^O rotational excitation temperature of 18 
K towards IRAS 13039, very close to the value we adopted. The 
approximate-equality T]^i„ ^ rex[C'^0(2 - 1)] was also recently 
found by Miettinen et al. (2011) for a sample of clumps asso- 
ciated with IRDCs. For the linear rotors SiO and DCN, we as- 
sumed that Tex is in the range [5 K, Eu/^b], whereas for CH3OH 
it was assumed that Tex e [10 K, 2Eu/3kB]- The adopted up- 
per Tex limits result from the turning point of the temperature- 
dependent part of A'(mol), i.e., from the derivative 



dA^(mol) 



-Zrot(rex)e 



0, 



(A.7) 



d^ex dFe;^ 

where it is assumed that Tex » Tbg (e.g.. |Hatchell et al. 19981 ). 
The upper Tgx limits thus derived lead to the lower limit to 
A^(mol). We note that the rex(CH30H) range we have used, 10- 
23.2 K, is comparable to the CH30H(2i: - 1^) rotational tempe- 
ratures of ~ 13.4 - 24.5 K derived by SSHIO towards a sample 
of massive clumps within IRDCs. 

We calculated the fractional abundances of the molecules by 
dividing the molecular column density by the H2 column den- 
sity: x(mol) - N{mol)/N(H2). For this purpose, the values of 
A^(H2) were derived from the LABOCA dust continuum map 
smoothed to the corresponding resolution of the line observa- 
tions. 



A.4. CO depletion factors 

Because the clumps under investigation (or at least most of them) 
are presumably cold, the CO molecules should efficiently freeze 
out onto dust grain surfaces. To estimate the amount of this CO 
depletion, we calculated the CO depletion factors, /d, following 
the analysis presented in the paper by Fontani et al. (2006). If 
x(CO)can is the "canonical" (undepleted) abundance of CO, and 
x(CO)obs is the observed CO abundance, /d is given by 



xjCO), 
jc(CO) 



(A.8) 



obs 



The "canonical" CO abundance at the galactocentric dis- 
tance Rqc was calculated using the relationship [Eq. (7) in 
IFontanietal. 20061 



x(CO)can = 9.5 X 10 



-5gl.l05-0.13«Gc[kpc] 



(A.9) 
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Table A.l. Spectroscopic properties of the observed transitions. 



Transition 


[cm J 


Eu/ks 
Li^J 


'1 c 

LiJ J 


A 
Livirizj 


B 

[iVlrlZJ 


C 

[iVlrlZJ 


SiO(5 - 4) 


2.3 X 10" 


31.3 


48.14651 




21711.96 




DCN(3 - 2) 


~ 10' 


20.9 


80.50709 




36 207.46 




"CO(2 - 1) 


8.9 X 10^ 


15.9 


0.02437 




55 101.01 




C"0(2- 1) 


9.5 X 10^ 


16.2 


0.02432 




56 179.99 




CH,OH(5_i5 -4_i4)-E 


9.5 X 10^ 


40.4 


3.88240 


127 484.0 


24 679.98 


23 769.70 


CH3 0H(5o,5-4o,4)-A+ 


1.1 X 10" 


34.8 


4.04297 




-II- 


-II- 


CH,OH(54,. -44,.)-A+/-'' 


3.5 X 10^ 


115.2 


1.45278 




- II- 





Notes. Column (2) gives the critical density of the transition, n^^t = Aa\ICu\, where Aui is the spontaneous decay rate and Ca\ is the collisional 
de-excitation rate. It was calculated at T = 15 K using the coUisional-rate data available in the Leiden Atomic and Molecular Database 
(LAMDA; |http : // www . strw . leidenuniv . nl/~moldata/ I dSchoier et al. 2005T l. For DCN, there is no molecular datafile available, so we 
used the collisional rate coefficient of HCN to estimate n^rit. Columns (3) and (4) give the upper-state energy and the product ^i\^S , where 
^ei is the permanent electric dipole moment, and 5 is the line strength. Columns (5)~(7) list the rotational constants of the molecules. The 
data were compiled from the Jet Propulsion Laboratory (JPL; IPickett et al. 1998t and CDMS JMiiller et al. 2005b databases.'"' A blend of 
CH30H(54,i - 44,o)-A+ and CH30H(54,2 - 44j)-A", which have the same spectroscopic parameters. 



Note that this relationship results from using the value = 8.5 
kpc, whereas our /?gc value was computed using - 8.4 kpc. 
Such a small discrepancy is however negligible. At Rcc - 8.5 
kpc, the above relationship gives the standard value 9.5 x 10"^ 
for the abundance of the main CO isotopologue in the solar 



neighbourhood ( Frerking et al. 1982| . At the galactocentric dis- 



tance of G304.74, x(CO)can - 11 X 10 To calculate the 
"ca nonical" C'^O abundan ce we employ the oxygen-isotopic x?i- 
tio ( I Wilson & Rood 1994| i 

= 58.8 X /?Gc[kpc] + 37.1 , (A.IO) 
and the ratio ["^0]/[''0] = 3.52 ( [Frerking et al. 1982| l: 



rCC'^Ol = x{CO\^n 

XK^ U;ea„ [180]/[170]X[160]/[180] 



x(CO)e.„ ^^^^^^ 



3.52 X (58.8 x/^GC [kpc] H- 37.1) 

For G304.74, we get x(C''0)ca„ ^ Jc(CO)can/1633 ^ 6.8 x 

10"^. The depletion factor /d is then calculated from /d - 

4C'^0)ea„/x(C''0)obs. 
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